Alloying-related trends from first principles: An application to the Ti-Al-X— N system 
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Tailoring and improving material properties by alloying is a long-known and used concept. Recent 
research has demonstrated the potential of ah initio calculations in understanding the material 
properties at the nanoscale. Here we present a systematic overview of alloying trends when early- 
transition metals (Y, Zr, Nb, Hf, Ta) are added in the Tii-^Al^N system, routinely used as a 
protective hard coating. The alloy lattice parameters tend to be larger than the corresponding 
linearised Vegard's estimation, with the largest deviation more than 2.5% obtained for Y0.5AI0.5N. 
The chemical strengthening is most pronounced for Ta and Nb, although also causing smallest elastic 
distortions of the lattice due to their atomic radii being comparable with Ti and Al. This is further 
supported by the analysis of the electronic density of states. Finally, mixing enthalpy as a measure 
of the driving force for decomposition into the stable constituents, is enhanced by adding Y, Zr and 
Nb, suggesting that the onset of spinodal decomposition will appear in these cases for lower thermal 
loads than for Hf and Ta alloyed Tii-^Al^N. 

PACS numbers: 61.66.Dk 68.60.Dv 71.15.Mb 71.20.Be 81.05.Je 
Keywords: Density Functional Theory; alloys; TiAIN; phase stability 
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I. INTRODUCTION 

Tii-^Al^N is nowadays a well-established material 
used as a protective hard coating due to its high hardness, 
relatively low coefficient of friction, and good oxidation 
and corrosion resistance [HQ- As such it has attracted 
significant attention both in basic as well as in applied 
research areas, resulting also in several quantum mechan- 
ical studies on phase stability, mechanical and electronic 
properties of the ternary Tii_ x Al x N system [§4ll|. 

One of the possible ways how to further tune the ma- 
terial properties is the concept of multicomponent alloy- 
ing. For example, Y, Zr, and Hf have been shown to 
improve the oxidation resistance of Tii-^Al^N (l2l - fF| 
and Cri-^AlxN coatings Qa|. Ta to increase the high 
temperature durability [lfj, Cr, Nb, and Ta to retard 
the decomposition process to higher thermal loads [l7i — 
[l9l | . The most recent successful approach is to combine 
theoretical studies with experimental work in order to 
gain deeper understan ding of the exp erimental observa- 
tions (e.g., Refs. [0, EMI EBia). Good examples 
demonstrating the ability of ab initio calculations for 
guiding the experiment by predicting chemistry-related 
trends are e.g. Refs. @, U, HI for various TM-A1-N 
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ternary systems. However, similar systematic theoreti- 
cal study showing the alloying-related trends for quater- 
nary Tii-jc-j/Al^XyN (or, more precisely, pseudo-ternary 
TiN-AlN-XN) systems, is missing. 

In this paper we demonstrate the potential of the first 
principle calculations for providing such systematic and 
exhaustive information for Tii-^ALjN alloyed with Y, 
Zr, Hf, Nb, and Ta throughout the whole quaternary 
phase field. It is clearly shown that with some care, the 
calculated results can serve as reliable "trend-givers" for 
the materials design. Lastly, it should be noted that 
in the paper we discuss only some related topics like the 
ground state properties, chemical strengthening, or onset 
of the spinodal decomposition; others, such as the final 
decomposition phase (precipitation of A1N in its stable 
wurtzitc (B4, ZnS prototype) structure or influence on 
oxidation resistance, remain for the future studies. 



II. COMPUTATIONAL APPROACH 

We used Vienna Ab initio Simulation Package [23L [24| 
employing the Density Functional Theory [25|, [2(| to per- 
form quantum mechanical calculations. The exchange 
and correlation effects were described with Generalised 
Gradient Approximation as parametrised by Wang and 
Perdew (27j . and implemented in projector augmented 
plane- wave pseudo-potentials [28| . We used the supercell 
approach to model random alloys, in particular 3x3x2 
(containing 36 atoms) supcrcclls were constructed for the 
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cubic rock-salt Bl structures. To do so, we optimised the 
short-range order parameters for pairs up to the fourth 
nearest neighbour distance, for triplets up to the third 
nearest neighbour distance, and for quadruplets up to the 
second nearest neighbour distance, following the Special 
Quasi-random Structures methodology J^. The same 
technique was applied also to construct 2x2x2 super- 
cells (32 atoms) for the wurtzite B4 structures. More 
details on these cells can be found in Ref. (l3j . 

The plane wave cut-off of 500 eV and a minimum of 
ps 3500 k-point • atoms guarantee the total energy accu- 
racy in the order of meV per atom. The stability of 
systems can be described using energy of formation, Ef, 
calculated as 



1 r 



E f = ^(Ti^^ALXyN) - - (1 - x - y)E(Ti hcp ) 



xE(Al icc ) + yE(Xt) + -E(N 2 : 



(1) 



Here, £'(Tii_. r _ y AL c Xj / N) is the total energy per atom of 
c-Tii^-yAUXyN, £(Ti hcp ), £(Al fcc ) and E(X*) are the 
total energies of Ti in hexagonal close-packed (hep, A3), 
Al in face-centered cubic (fee, Al), X=Y, Zr, and Hf in 
hep, and X=Nb, and Ta in body-centered cubic (bec, A2) 
structures, respectively. E(N 2 ) denotes the total energy 
of a nitrogen molecule. 



III. RESULTS AND DISCUSSION 

A. Bulk properties and phase stability of 
quaternary Ti— Al— X-N 

Fitting the energy-volume curve with the Birch- 
Murnaghan equation of state (EOS) [3(| yields the 
ground state properties: lattice parameter, a, total en- 
ergy, Eg-, and bulk modulus, B. The cubic lattice pa- 
rameters of five quaternary systems investigated here as 
a function of A1N mole fraction, x, and XN mole frac- 
tion, y, exhibit almost linear behaviour. Straight equally 
spaced contours demonstrate this in Fig.[TJi on the exam- 
ple of Tii-^-yAl^HfyN. The linearly interpolated lattice 
parameter according to the Vegard's empirical rule [3l| 
reads 



a v (x, y) = xclaw + ya X N + (1 - x - y)a Ta 



(2) 



The difference Aa = a — ay is shown in Fig. [TJd again 
for the Tii-z-yALjHfyN system. It follows that the 
linear Vegard's estimate is erroneous by as much as 
1.2% in the middle of the HfN-AIN tie-line. Similarly, 
differences of ss 1.5%, ~ 1.0%, and w 1.0% are ob- 
tained also for Tii-^-yAl^ZryN, Tii-^-yAl^NbyN, and 
Tii-^-yAl^Ta^N systems, respectively, with the maxi- 
mum appearing always in the middle of the A1N-XN tie- 
line (X=Zr, Nb, Ta). Tii-s-^AlaY^N exhibits an error 
of > 2.5% for Y0.5AI0.5N, reflecting the fact that Y is 
a significantly larger atom than Al. These deviations 



should be considered when using the Vegard's linearisa- 
tion as an estimate of a. The reason for the non-linear 
behaviour of a can be traced back to the gradual changes 
in the bonding (e.g., from ionic to metallic character), as 
demonstrated for Tii_ x AL;N in Ref. Q. 

In all five cases of Ti^- x -yAl x 'KyN, alloying an element 
X to Tii-^ALzN increases the lattice constant when the 
Al content on the metallic sublattice, x, is kept constant 
as well as for the constant Al-to-Ti ratio, x/(l — x — y). 
These trends are in a good agreement with previous 
experimental observations and theoretical studies 
Il4l . Il8l - [2H . The calculated lattice parameters are sum- 
marised in Table HI 

Bulk moduli, being the measure of volume compress- 
ibility (the elastic behaviour) of the system, are plot- 
ted in Fig. [5] In all cases, the variation is non-linear 
but smooth. The results can be classified into three 
different compositional dependencies, based on the va- 
lency of the alloying element X. For the group IVB el- 
ements, Zr and Hf, it seems that the main parame- 
ter controlling B is the TiN mole fraction, z = 1 — 
x — y (there is only a minor variation along the tie- 
lines with a fixed TiN mole fraction). For exam- 
ple, bulk modulus of Tii_ T _yAl x HfyN varies between 
« 240 GPa obtained for Al-excess Hfi.^Al^N, and f=s 
290 GPa for TiN H (see Fig. [5b). For the isova- 
lent alloy Tii-^-yAl^ZryN, a similar behaviour is ob- 
tained (Fig. [2jl). Situation is somewhat different for 
the Tii-^-yALNbyN and Tii-^-yALjTayN alloys (group 
VB), as NbN (B = 307 GPa) and TaN (B = 330 GPa) 
[32^ . respectively, yield the highest bulk moduli of the 
quaternary systems investigated (Figs. [2t and d). For 
these systems, the main parameter controlling B is the 
A1N mole fraction, x. Finally, yet a completely differ- 
ent elastic response is obtained for Tii-^-yALcYyN alloy 
(.Byn = 130 GPa), where the bulk modulus reaches val- 
ues well below 100 GPa for compositions in the centre 
of the compositional triangle (Fig. [5^). This is likely to 
be related to a strong instability of such material due to 
huge differences in the atomic size of individual species. 

A comparison of Figs.[T^ and[5Ji together with the data 
in Table |T] offers the possibility to design lattice matched 
superlattice materials (i.e. multilayer arrangements with 
layer thicknesses within the order of the lattice parame- 
ter) with variations in bulk modulus. For example, TiN 
and Hfo.3 AI0.7N have the same equilibrium lattice param- 
eters of ss 4.25 A, but about 40 GPa difference in B. An- 
other example is Tio.10Alo.3sZro.52N and TaN, both hav- 
ing lattice parameter a « 4.42 A, but the respective bulk 
moduli being 230 GPa and 330 GPa. Systems consisting 
of varying softer and harder components are known to be 
able to e.g., hinder dislocation motion, or stop or deflect 
propagation of cracks [H, HH . Here, an additional ben- 
efit stems from expected perfect interfaces (same crystal 
structures, lattice- matched) . 

It has been shown that for cubic structured Tii-^ALN 
thin films, the hardness reaches a maximum for the high- 
est Al contents for which the cubic structure is still main- 
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FIG. 1. a) Lattice parameter, a, and b) its deviation, Aa, from the Vegard's estimate for the quasi-ternary TiN-AlN-HfN 
system. 
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FIG. 2. Bulk modulus, B, as a function of the composition of pseudo-ternary alloys: (a) TizAbY^N, (b) TizAbZr^N , (c) 
TizALcNb^N, (d) TizALNbyN, and (e) TizAbTa y N. Countours are shown every 25 GPa for Ti 2 AbY^N (a), and every 5 GPa 
for the other four systems (b)-(e). 



tained [35||. The chemical compositions with hardness 
maximum are close to Tio.34Alo.66N for physical vapor- 
deposited coatings. These compositions provide also the 
highest driving force for decomposing towards their sta- 
ble constituents, cubic TiN and wurtzite AIN, via a spin- 
odal decomposition route. The calculated energies of for- 
mation, Ef , of the cubic and wurtzite phases were inter- 
polated with cubic polynomials in x and y, as imple- 
mented in the software package Mathematica. Equating 
£^. ub = EJ ur yields the cross-over, an estimation of the 
cubic and wurtzite single phase fields. Considering only 



small amounts of X (i.e., Y, Zr, Nb, Hf, or Ta) up to rj 0.1 
on the metallic sublattice, Nb and Ta do not change the 
maximum solubility of Al in the cubic phase consider- 
ably [H, lU, Zr and Hf decrease it slightly [H 0, [H, 
while Y decreases the Al solubility limit within the cubic 
structure drastically from x w 0.7 to x w 0.56 [l2l. I20I]. 
The lowering of the solubility limit in the case of Zr and 
Hf is a result of widening of the dual-phase region in the 
Zri_ K Al x N and Hfi-^Al^N phase diagrams, as discussed 
in Ref . Q . The lowered solubility limit of Y is likely to be 
connected with the large differences in the atomic sizes 
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FIG. 3. Energy of formation, Ef as a function of the va- 
lence electron concentration (VEC) for Tio.5Alo.5N, Ti-excess 
(Ti0.5Al0.444X0.056N), and Al-excess (Ti0.444Al0.5X0.056N) al- 
loys. 



(r A i = 1.432 A, r Ti = 1.448 A, and r Y = 1.80 A) destabil- 
ising the close packed rock-salt structure while favouring 
the more open wurtzite structure. 



B. Alloying effects on Tio.5Alo.5N 

In order to elucidate the role of individual ele- 
ments, we use c-Tio.5Alo.5N as a test composition. 
By replacing 1 Al or 1 Ti (out of 18 on the metal- 
lic sublattice) with 1 X atom, the alloy composi- 
tion is changed to Tio.sooAlo. 444X0. 056^ (Ti-excess) or 
Ti0.444Al0.500X0.056N (Al-excess) . 



1. Chemical strength and elastic distortions 

The relative phase stability change by alloying may 
be quantified by energy of formation expressing the en- 
ergy gain of forming an alloy of a given composition with 
respect to the individual elements. The results are plot- 
ted in Fig. [3] against the valence electron concentration 
(VEC) which is calculated as an average value of va- 
lence electrons per formula unit (f.u., 1 metal and 1 ni- 
trogen atom). Ti ([Ar]3d 2 4s 2 ), Al ([NcjSs^p 1 ), and N 
([Hc]2,s 2 2p 3 ) have 4, 3, and 5 valence electrons, respec- 
tively, hence VEC of Tio.5Alo.5N is i(|(4 + 3) + 5) = 8.5. 
These results indicate that Y, Ta, Nb, and Zr slightly 
increase the energy of formation while Hf works in the 
opposite direction. Nevertheless, these changes are not 
large (below 50meV/at.) and could be influenced by the 
actual supercell configuration. In addition, one should 
bear in mind that these are results of K calculations, 



and thus the mutual differences may change at higher 
temperatures due to, e.g. vibrational contribution to the 
Gibbs free energy. 

In order to address the strengthening effect of the al- 
loying elements, we first introduce cohesive energy calcu- 
lated as 

E c = E(Th_ x _ y Al x XyN) - i [(1 - x - y).E(Ti) 

+ xE(Al) + yE(X) + £(N)1 (3) 

where E(X) is the total energy of an isolated atom X. 
Based on this quantity, two energy terms are discussed: 
relative chemical strength, calculated as the cohesive en- 
ergy difference between relaxed Tio.5Alo.5N and (unre- 
laxcd) Tio.5Alo.5N+lX, where either one Ti or Al atom 
(out of 18 in the supercell) is replaced with one X atom. 
In the second step we allow for a full structural relax- 
ation. The total energy (or cohesive energy) decrease 
due to the relaxation (with fixed chemical composition 
of the alloy) corresponds to the elastic energy related to 
the local distortions caused by the foreign atom. These 
two contributions are shown in Figs. 3ji and b. 

The relative chemical strengthening (Fig. [4^,) demon- 
strates the detrimental influence of Y. For all the 
other elements the relative chemical strength of 
(Tio.5Alo.5)i-j,Xj,N with respect to Tio.5Alo.5N becomes 
more negative. The cohesive energy, as defined in Eq. |3l 
expresses the energy (per atom) needed to de-assemble 
the crystalline material into individual atoms. Since the 
clastic relaxation of the atom positions can only lower 
the total or cohesive energy, a negative/positive value 
of the relative chemical strengthening is an indicator for 
stronger/ weaker bonded atoms in the crystal. This is 
a purely electronic effect (a local change of the elec- 
tronic structure) and as such we interpreted it as elec- 
tronic strengthening. Consequently, the quaternary al- 
loys with Zr, Hf, Nb, and Ta are expected to be more 
resistant against plastic deformation, leading e.g., to a 
change in hardness of the alloys (although there are many 
other aspects influencing material's hardness). Indeed, 
a comparison with the experimental values taken from 
the literature, -J/tiAin ~ 33GPa [Hj], -ffnAiYN ~ 23GPa 
0, #TiAiZrN « 38 GPa 0, tfxiAiNbN « 37 GPa M, 
^TiAiHtN « 35 GPa 0, and # T iAiTaN « 37 GPa [ljj, 



supports this qualitative trend. It should be mentioned 
that the experimental value of -Htiaiyn corresponds to 
a wurtzite single phase structure, a phase change caused 
by the addition of w 10mol% YN to Tio.5Alo.5N. Since 
the wurtzite phase is in general softer than the cubic one, 
this significant decrease in hardness is likely to combine 
both effect, the phase change as well as the weakening of 
the chemical bonding due to Y. Finally, the linear guide- 
for-the eye suggest that the chemical strength increases 
with increasing VEC, at least in the investigated region. 

Figure 0)d shows how much elastic energy is intro- 
duced by local distortions into the material by alloy- 
ing. The elastic energy is plotted against the volume 



5 



(a) 



„ 6 



4 - 



CD 
LD 
Lfj 

> 

CD 



D) 
CD " 

v> 

TO 
O 



Al-excess 



Ti-excess 



Al-excess 



* 


Ti 0.5 AI 0.5 N 


< 


+Y 


■ 


+Zr 


► 


+Nb 


▲ 


+Hf 


T 


+Ta 



Al-e? 



-4 - 



8.4 



8.45 




(b) 




9.3 9.4 
volume [A 3 /at.] 



FIG. 4. (a) Chemical strengthening and (b) alloying- 
related elastic energy caused by the local distortions of 
Tio.5Alo.5N, Ti-excess (Ti0.5Al0.444X0.056N), and Al-excess 
(Ti0.444Al0.5X0.056N) alloys. 



per atom. The volume expansion in the first approxima- 
tion equals to 3e, where e is the homogeneous isotropic 
strain. Assuming the elastic constants are not signifi- 
cantly influenced by alloying, the elastic energy is pro- 
portional to e 2 . The quadratic fit gives qualitatively 
a good description as can be seen in Fig. 2b. Taking 
into account the metallic radii of individual metal atoms, 
tai = 1.432 A and tt; = 1.448 A, allows to explain why 
the smallest elastic deformations are caused by alloying 
with Nb (r Nb = 1.46 A) and Ta (r Ta = 1.46 A), fol- 



lowed by Hf (r Hf = 1.564 A) and Zr (r Zr = 1.59 A), and 
the largest distortions being caused by alloying with Y 
(ry = 1.80 A), as in this sequence the atomic radii in- 
crease. 

In conclusion, Ta and Nb are predicted to be the most 
dominant strengthening elements, which at the same 
time introduce also the smallest lattice distortions. 



2. Density of states 

The total density of states (DOS) for Tio.5Alo.5N and 
corresponding Al-excess alloys is shown in Fig. [5p,. The 
top of the valence band, as shown, can be divided into 
three distinct regions: (i) —10 to « — 6eV corresponding 
to hybridised sp 3 d 2 covalent bonding (as demonstrated 
by the overlap of ,s-, p— , and ci-projected DOS in Fig.[5ji), 
(ii) -6 to « — 2eV corresponding predominantly to hy- 
bridised N-p-TM-ci bonding, and (iii) w —2 to eV corre- 
sponding to metallic bonding formed by TM d electrons. 
The term "covalent" in the following analysis refers to the 
regions (i) and (ii) together (i.e., energy between —10 and 
ss -2cV). 

To quantify the effect of alloying elements, we intro- 
duce "centre of mass" (COM) [19|, [2(| of the above men- 
tioned energy bands, calculated as average energy of the 
band weighted by DOS: 



COM 



(4) 



where p{E) is DOS in the band between E n 
In assessing the alloying effects, the relative decrease 
of COM in particular energy range contributes towards 
lower total energy, and consequently leads to chemical 
strengthening as discussed in Section MI B II Addition- 
ally, each energy range (or band) can be characterised by 
a "band width" (BW) defined as: 



BW 



\ fE^P(E)dE 



(5) 



Such quantity reflects the width of the band weighted by 
the density of states distribution. Smaller value of BW 
corresponds to better overlap of the hybridised states, 
and hence to increased degree of hybridisation. 

COM of the metallic range of Tio.5Alo.5N is -0.684 eV. 
This value changes to -0.663, -0.688, -0.768, -0.662, 
and — 0.762 eV for Y, Zr, Nb, Hf, and Ta, respectively, 
as evaluated for the Al-excess compositions (see Fig.O)). 
Here, a significant strengthening is predicted by alloying 
of Nb and Ta, related to an increased concentration of 
d electrons forming the d-d bonds (t2 g symmetry) along 
(110) directions 0, [2(|. The same analysis of COM for 
the covalent region yields — 5.148 eV for Tio.5Alo.5N, and 
-5.034, -5.173, -5.192. -5.209, and -5.256 for Y, Zr, 
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FIG. 5. (a) Density of states for pure, Y-, Zr-, Hf-, Nb-, 
and Ta-containing Tio.5Alo.5N for the Al-excess compositions 
(y = 0.056). The s-, p-, and rf-projected density of states 
of Tio.5Alo.5N is shown on the very top. Ef denotes the 
Fermi energy, (b) Calculated corresponding centres of mass 
for metallic and convalent regions, and the band width of the 
sp 3 d 2 hybridised region (—10 and ~ — 2eV). 



Nb, Hf, and Ta, respectively, in the Al-excess configura- 
tions (Fig. [SJd). These numbers suggest that Y consid- 
erably weakens the covalent interaction, while the other 
four elements (Zr, Nb, Hf, and Ta) all strengthen the co- 
valent bonding. Finally, a consistent picture is obtained 
also from BW of the hybridised region: Y significantly 
broadens the electronic states distribution while Ta con- 
fines the states resulting in an increased degree of the 
sp 3 d 2 hybridisation. 

In summary, the most significant strengthening is ob- 
tained by alloying with Ta, followed by Nb and Hf. Still 
a small strengthening is obtained by alloying with Zr, 
while Y softens the material. These trends are in line 
with the chemical strength considerations based on the 
total energy changes, as discussed in Section UlIB II 

The density of states at the Fermi level is an indica- 
tor of the (relative) alloy stability. It reaches a value of 
0.21st.-cV _1 -at. _1 for Tio.5Alo.5N. Zr keeps it unaltered 
for the Al-excess composition, while Hf, Nb, Ta, and Y 
increase it to 0.24, 0.25, 0.25, and 0.29 st. • eV" 1 • at." 1 , 
suggesting a gradual destabilisation of Tii-^Al^N by al- 
loying. This is in qualitative agreement with the trends 
in energy of formation as shown in Fig. [3] 



3. Decomposition of the unstable alloys 

The mixing enthalpy, i? m i x , is calculated as 

#mix = S(Ti 1 _ K _ y Al x X a N) - [(1 - x - y)E(TiN) 

+ xE(AlN) + yE(XN)] , (6) 

where E(XN) is the total energy (or energy of for- 
mation, or cohesive energy) of binary cubic XN, and 
E(Tii_ x _yXl x Xy~N) is the corresponding energy of the 
quaternary alloy. The mixing enthalpy expresses the en- 
ergy gain (when negative) of forming the alloy with re- 
spect to the binary nitrides, hence it quantifies whether 
the alloy is stable (i? m i x < 0) or unstable (H m - lx > 0) 
with respect to decomposition into the cubic binary ni- 
trides. 

The results obtained are plotted in Fig.[S]for the whole 
compositional range 0<x<l,0<j/<l, and 
< z < 1, x+y+z = 1, for all five quasi-ternary systems. 
The driving force for decomposition is the largest for 
Tii_. c _.yAL s Y. y N near to the YN-A1N tie-line (Fig. [5Ji). 
The lattice difference between YN and A1N is about 19%, 
a value larger than what is "allowed" for solid solutions 
by Hume-Rothery rules . H m i x is positive in the whole 
compositional range except for alloys close to the quasi- 
binary TiN-NbN and TiN-TaN tie-line (Figs. Efc and c). 
This suggests that the Tii_ y Nb y N and Tii_ y Ta y N al- 
loys are stable, a result previously shown in literature 
pH [37l HU . The calculated mixing enthalpies are sum- 
marised in Table UJ 

Adding a small amount (up to y = 0.1) of Y, Zr, or Hf 
to Tio.5Alo.5N while either keeping the Al-to-Ti ratio, or 
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keeping the Al amount x = 0.5 constant, the driving force 
for decomposition H m i X increases. Tii-^Al^N has been 
shown to decompose spinodally before the stable wurtzite 
A1N precipitates appear [39[ , and also the quaternary al- 
loys are expected to do so. Therefore, the increased force 
for decomposition is interpreted as increased force for 
the isostructural decomposition, thus shifting the onset 
of spinodal decomposition to lower temperatures (lower 
thermal loads) . Consequently, age- hardening of these al- 
loys is predicted to occur at lower temperatures than 
for Tio.5Alo.5N. The peak broadening in X-ray diffrac- 
tograms suggesting the spinodal decomposition is indeed 
observed for lower annealing temperatures when Zr fl3j 
is added . On the contrary, addition of Nb and Ta lowers 
the mixing enthalpy thus the onset of the spinodal de- 
composition is predicted to be shifted to slightly higher 
thermal loads as compared with Tio.5Alo.5N, as in excel- 
lent agreement with experiments, see e.g. Ref. (l9l |. 



viation from Vegard's-like linear relationship, being the 
largest with w 2.5% for Y0.5AI0.5N. The additional com- 
positional degree of freedom of the quaternary alloys (as 
compared with Tii-^Al^N) is suggested to allow for spe- 
cially designed lattice-matched multilayers with alternat- 
ing soft and hard components. 

Analysis of the chemical strength and local elastic dis- 
tortions showed that Ta and Nb are the most promising 
strengtheners, closely followed by Zr and Hf. These con- 
clusions agree well with the published experimental re- 
sults of the increased hardness of Tii-^Al^N alloyed with 
Zr, Nb, Hf, or Ta. Finally, the analysis of the mixing 
enthalpy as a measure of the driving force for decompo- 
sition suggests that addition of Y, Zr and Hf leads to its 
increase, thus earlier onset of the spinodal decomposition 
and the related age-hardening process. 
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Appendix A: Properties of the quasi-ternary alloys 
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TABLE I. Calculated lattice parameter, a, bulk modulus, B, and mixing enthalpy, //mix, of quasi-ternary TiN-AlN-XN systems for X =Y, Zr, Nb, Hf, and Ta, as a 
function of TiN (z), A1N (x) and XN (y) mole fractions. 
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